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Abstract 

Understanding the electronic properties of quasicrystals, in particular the dependence of these 
properties on dimension, is among the interesting open problems in the field of quasicrystals. We 
investigate an off-diagonal tight-binding hamiltonian on the separable square and cubic Fibonacci 
quasicrystals. We use the well-studied singular-continuous energy spectrum of the 1-dimensional 
Fibonacci quasicrystal to obtain exact results regarding the transitions between different spectral 
behaviors of the square and cubic quasicrystals. We use analytical results for the addition of Id 
spectra to obtain bounds on the range in which the higher-dimensional spectra contain an abso- 
lutely continuous component. We also perform a direct numerical study of the spectra, obtaining 
good results for the square Fibonacci quasicrystal, and rough estimates for the cubic Fibonacci 
quasicrystal. 

1 Background and Motivation 

As we celebrate the Silver Jubilee of the 1982 discovery of quasicrystals [1], and highlight the achieve- 
ments of the past two and a half decades of research on quasicrystals, we are reminded that there still 
remains a disturbing gap in our understanding of their electronic properties. Among the open questions 
is a lack of understanding of the dependence of electronic properties — such as the nature of electronic 
wave functions, their energy spectra, and the nature of electronic transport — on the dimension of the 
quasicrystal. In an attempt to bridge some of this gap, we [2,3] have been studying the spectrum and 
electronic wave functions of an off-diagonal tight-binding hamiltonian on the separable n-dimensional 
Fibonacci quasicrystal^] [4]. The advantage of using such separable models, despite the fact that they 
do not occur in nature, is the ability to obtain exact results in one, two, and three dimensions, and 
compare them directly to each other. Here we focus on the energy spectra of the 2-dimensional (2d) 
and 3-dimensional (3d) Fibonacci quasicrystals to obtain a quantitative understanding of the nature 
of the transitions between different spectral behaviors in these crystals, as their dimension increases 
from 1 up to 3. In particular, we consider the transitions between different regimes in the spectrum, 
taking into account the existence of a regime in which the spectrum contains both singular continuous 



1 The reader is referred to Refs. [5] and [6] for precise definitions of the terms 'crystal' and 'quasicrystal'. 
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and absolutely continuous components. These different behaviors of the higher-dimensional spectra 
are expected to reflect on the physical extent of the electronic wave functions, as well as on the dy- 
namics of electronic wave packets, and are therefore of great importance in unraveling the electronic 
properties of quasicrystals in general. 

Recall [2] that the off-diagonal tight-binding model assumes equal on-site energies (taken to be 
zero), and hopping that is restricted along tile edges, with amplitude 1 for long (L) edges and T for 
short (S) edges, where we take T > 1. The Schrodinger equation for the square Fibonacci quasicrystal 
in 2d (with obvious extensions to higher dimensions) is then given by 

T n+1 V(n + 1, m) + T n ^{n - 1, m) + T m+1 V(n, m + 1) + T m ^(n, m - 1) = EV(n, m), (1) 

where ^(n,m) is the value of a 2d eigenfunction on a vertex labeled by the two integers n and m, 
and E is the corresponding eigenvalue. The hopping amplitudes Tj are equal to 1 or T according to 
the Fibonacci sequence {!)} = {1, T, 1, 1, T, 1, T, 1, 1, T, 1, 1, T, 1, T, 1, 1, T, 1, T, . . .}. By prohibiting 
diagonal hopping, the resulting high-dimensional eigenvalue problem is ensured to be separable. This 
allows one to use the known solutions for the Id problem [7-14] in order to construct the solutions 
in two and higher dimensions (as was done for similar models in the past [15-21]). Two-dimensional 
eigenfunctions can therefore be expressed as Cartesian products of the Id eigenfunctions [3], and the 
corresponding 2d eigenvalues are given by pairwise sums of the Id eigenvalues. 

The Id spectrum for the N th order Fibonacci approximant is composed of Fjy bands, where Fjy = 
Fiv-i + Fn-2 is the N th Fibonacci number, starting with Fq = F\ = 1. The edges of each such band 
correspond to either periodic or antiperiodic boundary conditions. Hence, by direct diagonalization 
of the two corresponding hamiltonians for a single approximant we obtain the edges of the energy 
intervals in the spectrum. The 2d and 3d spectra are then calculated as the Minkowski sums of two or 
three Id spectra, where the Minkowski sum of two sets A and B is the result of adding every element 
of A to every element of B, i.e. the set 

A + B = {x + y | x e A, y € B} . (2) 

Although the spectrum of the Id Fibonacci model, for any choice of T ^ 1, is a totally disconnected 
set with zero bandwidth and an infinite number of bands, the higher-dimensional spectra exhibit 
different behavior for different values of the relative hopping parameter T, including spectra that 
contain continuous intervals and have a finite measure [2]. A similar situation arises in the case of 
the well-known ternary Cantor set [22], which is constructed iteratively by starting with the closed 
interval [0, 1], and at each iteration removing the open middle thirds of all remaining closed intervals 
from the previous iteration. The first few approximants that are obtained in this way are Cq = [0, 1], 
Ci = [0, l/3]U[2/3, 1], andC 2 = [0, 1/9] U [2/9, 1/3] U [2/3, 7/9] U [8/9, 1], so that after N such iterations 
one is left with an approximant set Cn consisting of 2 N closed intervals, each of which has a measure 
1/3^, and therefore the total measure of the set is (2/3) N . The ternary Cantor set itself Coo, defined 
as the limit — > oo of this sequence of sets, contains uncountably-many points yet no interval, it is 
totally disconnected, and its total measure is zero. By simple inspection, one finds that for any finite 
order Cantor approximant CV, the Minkowski sum Cn + Cn is the entire interval [0,2]. One can 
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show that this also holds in the limit N — * oo, namely that Coo + Coo = [0, 2]. Thus, even though 
contains no interval, its sum with itself covers the whole interval from to 2. 

For a given dimension n, we identify a sequence of values 1 < lf d) < T 2 (nd) < T 3 (nd) < T 4 (nd) 
corresponding to the following transitions in the spectrum: 

rp( nd ). rpj^ va j ue Q f p below which all bands in the n-dimensional spectrum are of positive, finite 
measure. For T > T^ n< ^ there is at least a finite number of zero measure bands in the spectrum. 

rp( nd ). rpj^ va j ue Q £ p a )3 0ve w hich the number of bands in the n-dimensional spectrum is infinite. An 
infinite number of bands in a spectrum of finite bandwidth necessarily implies that infinitely 
many bands are of zero measure, thus T^"^ > T^ nd ^ . 

rp( nd ). va j ue Q £ p a bove which all bands in the spectrum are of zero measured 

rp( nd ). rpk e va j ue Q f p ^Qyg w bich the total bandwidth of the spectrum is zero. 

We use two different approaches to study the behavior of the spectrum. In Sec. [2] we use analytical 
results derived for the addition of generalized Cantor sets to obtain an upper bound on the transition 
between absolutely continuous and singular continuous spectra. In Sec. [3] we use direct numerical 
calculation of the 2d and 3d spectra of Fibonacci approximants of finite order to extrapolate for the 
behavior in the quasiperiodic limit. In an earlier paper [3] we studied only two of the transitions, 
T^® and T^ nd . To find T^ nd ^ we used a naive method based on the results of Ashraff et al. [21] 
for the diagonal tight-binding hamiltonian. The current results include a correction to our previous 
calculation. In Sec. 0] we summarize the results, and discuss their expected relation to the nature of 
eigenfunctions and to quantum dynamics, indicating directions for future work. 



2 Analytical bounds for the appearance of continuous intervals in 
the spectrum 

2.1 Addition of generalized Cantor sets - Known results 

A generalized Cantor set is obtained just like the ternary Cantor set except that the open intervals 
removed at each iteration are not necessarily the middle thirds of the remaining closed intervals. For 
each interval removed from the set, one defines a left (right) ratio of dissection as the ratio between 
the length of the left (right) remaining interval and the length of the original one. Sets for which the 
left and right ratios are the same are called central Cantor sets. In general, the ratios of dissection 
may vary between the left and right resulting intervals, between different iterations of the process, 
and between different intervals at the same step. The ternary Cantor set is a central Cantor set with 
a constant ratio of dissection of 1/3. 

We are interested in conditions for the appearance of intervals in the Minkowski sum of n general- 
ized Cantor sets. For central Cantor sets with a constant ratio of dissection a, one can show that the 

2 Note that the absence of intervals in the spectrum above does not necessarily correspond to zero total band- 

width. It is in fact possible to use the Cantor set generation process to obtain a totally disconnected set with a finite 
measure. For example, if at the N th iteration of the generation process the middle 1/3^ part is removed from each of 
the remaining intervals, one ends with a totally disconnected set whose measure is 

limnr=i (1 - V3 fe ) ^ 0.5601. 
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condition for the sum to be an interval is 



n > 1 or a > . (3) 

1-a ~ ~ n + 1 

Thus, the ternary Cantor set exactly has the critical value of o = 1/3 for which a sum of re = 2 central 
Cantor sets is an interval. Cabrelli et al. [23] found, more generally, a sufficient condition for the 
existence of an interval in the sum of n generalized Cantor sets, all of which can be constructed with 
a lower bound a on their ratios of dissection, which is given by 

(»" 1)7^13 +T^->1- (4) 
(1 — a) 6 1 — a 

2.2 Applying Cantor set results to the Fibonacci spectra 

Before using the results quoted above to analyze the Fibonacci spectra, we should note that there 
exist two important differences between the energy spectra Sjy of the N th order approximants of the 
Id Fibonacci quasicrystal, and finite approximants Cn of generalized Cantor sets. The spectrum 
Sn consists of Fn rather than 2 N energy intervals, and is not contained in the spectrum SV-i of 
the approximant of order N — 1. One should therefore take care in defining the spectrum S of the 
Fibonacci quasicrystal itself as the set of limit values for sequences of energies taken from consecutive 
spectra Sn of finite order approximants 

S={E= lim E n | E n € S n \ . (5) 

The fact that the number of bands in Sn is Fn rather than 2^ implies that the spectra cannot 
be constructed by the iterative process described above for generalized Cantor sets, and hence that 
the ratios of dissection cannot be defined. However, the spectrum Sn of a finite approximant can be 
padded with additional intervals which can be chosen in a manner that will not disturb the calculation, 
and will increase the number of intervals to 2^, as in the Cantor approximant. This allows to calculate 
backwards and define effective ratios of dissection. The additional intervals can be added on either, 
or both, ends of the spectrum. Thus, the effective ratios of dissection are not uniquely determined. 

We have tried using Eq. (j3J) to find a sufficient condition for the higher-dimensional spectra to 

contain an interval. This would provide a lower bound on T^ nd ^ — a value of T below which the 

condition is satisfied and the re-dimensional spectrum necessarily contains an interval. Unfortunately, 

as one studies the effective ratios of dissection defined for the Id spectrum it turns out that regardless 

of the way in which the approximant spectra are embedded in Cantor approximants, the ratios of 

dissection are not bounded away from zero, even for small values of T, as shown in Fig. Ufa). Hence, 

at this point we do not know how to use the condition of Cabrelli et al. to obtain a lower bound on 

rp(nd) 

1 3 

Nevertheless, by studying the maximal effective ratio of dissection we can obtain an upper bound 
for the value of Tg^ above which the higher-dimensional spectra do not contain an interval. Fig. H{b) 
shows the effective maximal ratio for approximants of order N = 5 (F/v = 8) to N = 14 (Fn = 610). 
It is evident that the maximal ratio of dissection rapidly converges as a function of the order of the 
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approximant, with almost no difference between the the curves for N = 9 and above. It is also of 
interest to note that the maximal ratio is independent of the way in which the approximant spectrum 
is embedded in a Cantor approximant. Values of T for which the maximal ratio of dissection fails to 
satisfy Eq. ([3]) imply that there is no portion of the Id spectrum which can lead to the existence of an 
interval in the higher-dimensional spectra. The maximal ratio of dissection becomes 1/3 at T ~ 3.15 
and 1/4 at T ~ 4.2. Thus, we expect to see the vanishing of intervals in the spectrum at a value of T 
below these upper bounds for 2d and 3d respectively. 



3 Direct study of the 2d and 3d spectra 

We now turn to the direct study of the higher-dimensional spectra. This is done by explicitly cal- 
culating the spectra for approximants of finite order. Each pair or triplet of energy bands in the 
Id spectrum is summed to yield a single band in the 2d or 3d spectrum, respectively. A set of F/v 
bands in the Id spectrum generates (F/v + l)F/v/2 bands in the corresponding 2d spectrum, and 
(F/v + 2) (F/v + l)F/v/6 bands in the 3d spectrum, with possible overlaps that decrease as T increases. 
Overlapping bands are merged into single energy intervals to obtain the actual structure of the higher- 
dimensional spectra. Note that we shall use the term 'bands' to refer to the continuous energy intervals 
in the spectra, even though strictly speaking they may be composed of different bands with overlapping 
energies. 

3.1 Measuring the smallest band in the spectrum to find 

For T > there is at least one zero-measure band in the spectrum. We therefore measure the 

smallest band B m i n and ask whether it vanishes in the limit of N — > oo. For T < T^ nd ^ the length of 

(nd) 

the smallest band is independent of the order N of the approximant. For T > it can be described 
by a power law B m i n cx F N T ^ with some positive exponent, a n (T). We locate by finding the 

value of T for which a n vanishes. Fig.[^b) clearly shows that 1.6 < < 1.8, and Fig.[2Jd) indicates 
that 2 < T^ d ^ < 2.6. Within these bounds, the width of the smallest band oscillates between the two 
different limiting behaviors. 

As T increases and the overlap of bands vanishes, the smallest band in the n-dimensional spectrum 
is expected to be n times the smallest band of the Id spectrum. Hence for high values of T the 
exponents a n (T) should be independent of the dimension, because the multiplicative factor of n only 
adds a constant term in the semi-logarithmic scale. Fig. [3] shows the extracted exponents a n (T), 
indicating that they indeed coincide for all values of T above T^ nd ^ . 

3.2 Counting the number of bands to find 

Next we count the number of bands #B in the spectrum and ask whether it tends to infinity or 

3 (T) 

remains finite as N increases. Again, we express this number as a power law of the form j^B oc F N n , 
expecting (3 n {T) to vanish for T < ■ For the Id Fibonacci quasicrystal #B\ = F/v oc t n , where 
r is the golden mean. In higher dimensions, as the overlap between bands decreases with increasing 
T, we expect the number of bands to tend to its maximal value, which is approximately (#i?i) 2 /2 in 
2d, and approximately (#-£>i) 3 /6 in 3d. Thus the exponents /3 n {T) should tend to nlogr as T — > oo. 
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The dashed horizontal line in Fig. [5] indicates the expected limit value for the 2c? model which indeed 
tends to it. For the 3d model the limit is only obtained at significantly higher values of T, indicating 
that the overlap of bands plays a significant role in the structure of the spectrum even at relatively 
high values of T. The continuous variation of (^(T) allows us to use smooth extrapolation and find 
rp(?d) ^ i qq^ whereas in 3d we can only conclude that 2.0 < < 2.6. Combining the fact that 

rp(nd) y rp(nd) ^ Q results for the exponents (3 n (T) as shown in Fig. [5l we find that at least in 2d 
and 3d, = T^ 1 , and hence that there 

only a finite number of zero-measure bands. 



and 3d, = T^ 1 , and hence that there is no intermediate regime in which the spectrum contains 



3.3 Measuring the largest band in the spectrum to find T^ nd ^ 

(rid) 

For T > Tg all bands in the spectrum have zero measure. We therefore look at the width of the 
largest band in the spectrum and ask whether it vanishes as N — > oo. However, since the maximal 
energy in the spectrum is approximately n(l + T), for small values of T the overlap of bands leads to 
an increase in the width of largest band as a function of T. To avoid this we normalize the results, 
dividing by the maximal energy in the spectrum. Thus, for T > T^ nd \ we express the normalized 
largest band as a power law B max oc F n ln ^ T \ Figs. [6^b) and [7] clearly indicate that T^ 2< ® ~ 2, but in 
3d oscillatory behavior dominates a large range of values for T, and we cannot determine the transition 
without extending the analysis to higher order approximants. However, from Fig. [6^d) we can infer 
that the transition occurs at some value of T below 5, for which we obtained analytically a stricter 
upper bound of Tg ^ < 4.2 as shown in Fig. QJb). 

As for B m i n , at large values of T, B max is also expected to be n times the largest band of the Id 
spectrum, and hence the exponents should be independent of dimension. The fact that this does not 
occur indicates, once again, that the overlap of bands is still significant for values of T as large as 6. 

3.4 Calculating the total measure of the spectrum to find 

To find T 4 (nd) we measure the total bandwidths of the spectra as N increases, normalizing by 1 + T, 

p / r p\ _____ 

and looking for a power law decay of the normalized bandwidth W oc F N . Fig. E(c) shows a 
decrease in the normalized total measure of the spectrum as a function of T in 2d, but Fig. E^d) shows 
the total measure in 3d to be almost independent of iV for any given value of T. Thus, although the 3d 
spectrum consists only of zero measure bands for values of T above 5, its total measure remains finite 
over the entire range of T values studied. The exponents S n (T) are shown in Fig. [H The transition 
to zero total bandwidth in 2d occurs at ~ 2.6. In 3d we can only say that T 4 (M) > 6. 



4 Summary and future work 

The results of Sections [2] and [3] are summarized as follows 





rp(nd) rp{nd) 
1 1 — 1 2 


rp(nd) 


Upper bound for T^ nd ^ 


rp(nd) 
J 4 


2d 


~ 1.66 


~ 2 


3.15 


~ 2.6 


3d 


2.0 - 2.6 


< 5 


4.2 


> 6 



(6) 



6 



The transitions between different regimes in the spectrum are expected to reflect on the physical 
properties of the Fibonacci quasicrystals, on the nature of eigenfunctions and on the dynamics of 
electronic wave packets. For values of T above the transition T^ 1 ^ the higher-dimensional spectra 
are similar to the Id spectrum in being totally disconnected, singular continuous sets, and hence 
the eigenfunctions are expected to be critical, and wave packets are expected to display sub-ballistic 
dynamics. Note that the last transition T^ nd ^ is of no consequence for this matter because the spectrum 
is purely singular continuous both above and below this value. For values of T below the lowest 
transition point T^ nd ^ = T^ d \ where the spectra are absolutely continuous we expect to find extended 
eigenfunctions, and wave packets are expected to display ballistic dynamics. For the intermediate range 
between these transitions the spectra contain both singular continuous and absolutely continuous 
parts, and therefore we expect to find mixed ballistic and sub-ballistic dynamics, and some of the 
wave functions to be extended. 

We intend to complement these studies by simulating the dynamics of electronic wave functions to 
find whether transition between ballistic and sub-ballistic dynamics occur at the points found here. We 
also intend to use the degeneracy of wave functions in the 2d Fibonacci quasicrystal (as hypothesized 
in Ref. [2]) to construct maximally extended wave functions, again, we expect to find some qualitative 
change in the nature of these wave functions near the transition points indicated above. 
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Figure 1: (a) The minimal effective ratio of dissection calculated for the Id spectra of approximants 
of order 12-14 for values of T up to 1.5. The sharp drop near T = 1 means that no value of T satisfies 
the sufficient condition for obtaining an interval in the higher-dimensional spectra, (b) The maximal 
effective ratio of dissection calculated for the Id spectra of approximants of order 5-14 for values of T 
up to 6. The horizontal dotted lines are drawn at 1/3 and 1/4 to indicate the upper bounds for the 
value of T at which no intervals are to appear in the 2d and 3d spectra respectively. 
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Figure 2: The length of the smallest band B m i n in the spectrum of the 2d (top) and the 3d (bottom) 
Fibonacci quasicrystals. The length of the smallest band is plotted on the left as a function of T for 
different approximants, and on the right as a function of N for different values of T. The linear slopes 
in the semi- logarithmic plots as a function of N indicate a power law behavior, B m i n oc T - Na «( T ), 
The exponents a n (T) are plotted in Fig. O 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

T 



Figure 3: The exponents a n (T) extracted from Fig. O The curves for Id, 2d, and 3d all coincide for 



values of T above the transition at T n 



(nd) 
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Figure 4: The number of bands #B in the spectrum of the 2d (top) and the 3d (bottom) Fibonacci 
quasicrystals. The number of bands is plotted on the left as a function of T for different approximants, 
and on the right as a function of N for different values of T. The linear slopes in the semi-logarithmic 
plots as a function of N indicate a power law behavior, j^B oc t n ^ t ' . The exponents P n (T) are 
plotted in Fig. [5j 




Figure 5: The exponents /3 n (T), extracted from Fig. [H The horizontal dashed line indicates the 
expected asymptotic value of 2 log r ~ 0.418 for the 2d quasicrystal. 
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Figure 6: The normalized length B max of the largest band in the spectrum of the 2d (top) and the 
3d (bottom) Fibonacci quasicrystals. B max is plotted on the left as a function of T for different 
approximants, and on the right as a function of N for different values of T. The linear slopes in 
the semi-logarithmic plots as a function of N indicate a power law behavior, B max oc r~ N ^ n<yT \ The 
exponent 72 (T) is plotted in Fig. [7J 



□ 1d 

x 2d 



Figure 7: The exponent 72(7") extracted from Fig. O and compared with 71 (T). The asymptotic 
behavior in which all three curves are expected to coincide is not observed for the values of T shown. 
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Figure 8: The normalized total bandwidth W of the spectrum of the 2d (top) and the 3d (bottom) 
Fibonacci quasicrystals. W is plotted on the left as a function of T for different approximants, and on 
the right as a function of N for different values of T. The linear slopes in the semi-logarithmic plots 
as a function of N indicate a power law behavior, W oc T ~ NS «( T ) . The exponents 5 n (T) are plotted in 
Fig. M 



□ 1d 
x 2d 
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Figure 9: The exponents 6 n (T), extracted from Fig. El In 3d it is evident that the transition to zero 
bandwidth does not occur within the studied range of T values. 
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